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In this work, we have assumed the generahzed Vaidya solution in Lovelock theory of gravity 
in (n + 2)-dimensions. It has been shown that Gauss-Bonnet gravity, dimensionally continued 



Lovelock gravity and pure Lovelock gravity can be constructed by suitable choice of parameters. 
We have investigated the occurrence of singularities formed by the gravitational collapse in above 
three particular forms of Lovelock theory of gravity. The dependence of the nature of singularity 
on the existence of radial null geodesic for Vaidya space-time has been specially considered. In all 
the three models, we have shown that the nature of singularities (naked singularity or black hole) 
■ completely depend on the parameters. Choices of various parameters are shown in tabular form. In 

\^ \ Gauss-Bonnet gravity theory, it can be concluded that the possibility of naked singularity increases 

00 . with increase in dimensions. In dimensionally continued Lovelock gravity, the naked singularity 

I is possible for odd dimensions for several values of parameters. In pure Lovelock gravity, only 

black hole forms due to the gravitational collapse for any values of parameters. It has been shown 
that when accretion is taking place on a collapsing object, it is highly unlikely to get a black 
hole. Finally on considering the phantom era in the expanding universe it is observed that there 
is no possibility of formation of a black hole if we are in the Gauss-Bonnet gravity considering the 
accreting procedure upon a collapsing object. 



1 Introduction 

Gravitational collapse is one of the most important problem in classical general relativity. The study of grav- 
itational collapse was started by Oppenheimer and Snyder (1939). They studied collapse of dust with a static 
Schwarzschild exterior while interior space-time is represented by Friedman like solution. One would like to 
know whether, and under what initial conditions, gravitational collapse results in black hole (BH) formation. 
In particular, one would like to know if there are physical collapse solutions that lead to naked singularities 
(NS). In last few years, there have been extensive studies on gravitational collapse in order to investigate the 
nature of the singularities. The study of gravitational collapse of spherically symmetric space-times led to many 
examples (Eardly et al, 1979; Christodoulou, 1984; Neuman, 1986; Waugh et al, 1986, Dwivedi et al, 1989) 
of NSs. There is no general theory of the nature and visibility of singularities. There do exist a number of 
exact solutions of the Einstein equation which admit, depending upon the initial data, BHs or NSs (Joshi, 1993, 
2000; Clarke, 1993; Wald, 1997; Jhingan, 1999; Singh, 1999). In particular the Vaidya solution (Vaidya, 1951) 
(contains outgoing radiation) is extensively used to show that the end state of collapse for regular initial data 
results in a NS. 

Harko et al (2000) have studied the gravitational collapse of strange matter and analyzed the condition for 
formation of a NS in the spherically symmetric Vaidya space-time. It has been shown that depending on the 
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initial distribution of density and velocity and on the constitutive nature of the collapsing matter, either a BH or 
a NS is formed. Santos and collaborators (Santos, 1984, 1985; de Oliveira et al, 1985, 1986, 1987, 1988) included 
the dissipation in the source by allowing radial heat flow (while the body undergoes radiating collapse). Ghosh 
and Deskar (2000, 2003) have considered collapse of a radiating star with a plane symmetric boundary and 
have concluded with some general remarks. Ghosh et al (2002) have discussed the study of collapse of radiating 
star in Vaidya space-time in {n + 2)-dimensions. Wang et al (1999) has generalized the Vaidya solution which 
include most of the known solutions to the Einstein equation such as anti-dc-Sittcr charged Vaidya solution. 
Husain solution has been used to study the formation of a BH with short hair (Brown et al, 1997) and can be 
considered as a generalization of Vaidya solution (Wang et al, 1999). Recently, Patil et al (2005, 2006) have 
studied the gravitational collapse of the Husain solution in four and five dimensional space-times and Debnath 
et al (2008) have studied the gravitational collapse of the Husain solution in (n + 2)-dimensional space-times 
for various types of equation of state. 

The Lovelock gravity theory (Lovelock, 1971) is a generalization of Einstein gravity theory in higher dimen- 
sional space-times, but the Lovelock theory is not a higher derivative gravity theory. The Lagrangian of the 
Lovelock gravity consists of the dimensionally extended Euler densities 



i=0 

where p < [{n + l)/2] (where, [x] denotes the integer part of the number arbitrary constants with 

dimension of [length]'^^^'^ , {n+2) is the space-time dimension and £j are the Euler densities of a (2i)-dimensional 
manifold (Cai et al, 2008) 



Here, the generalized delta function is totally antisymmetric in both sets of indices. Cq is set to one, there- 
fore the constant cq is just the cosmological constant. Ci gives us the usual scalar curvature term. If we set 
ci = 1, £2 gives just the Gauss-Bonnet (GB)term. So the Einstein-Gauss-Bonnet(EGB) gravity is a special 
case of Lovelock's theory of gravitation, whose Lagrangian just contains the first three terms in (1). 

Static spherically symmetric BH solutions can be found in Lovelock theory in the sense that a metric func- 
tion is determined by solving for a real root of a polynomial equation (Myers et al, 1988; Wheeler, 1986). 
More recently, static, non-spherically symmetric BH solutions have been also found in the Lovelock gravity 
(Cai, 2004). Spherically symmetric BH solutions in the GB gravity have been found and discussed in (Myers 
et al, 1988; Whiler, 1986; Boulware et al, 1985 ) and rotating GB BHs have been discussed in (Kim et al, 
2008). Some exact solutions for a Vaidya-like solution in the EGB gravity have been found in (Kobayashi, 2005; 
Maeda, 2006; Dominguez et al, 2006; Ghosh et al, 2008). Gravitational collapse and study of NS formation in di- 
mensionally continued Lovelock gravity have been investigated in ref.(Nozawa et al, 2006; Ilha et al, 1997, 1999). 

In this paper, we are mainly studying the nature of singularities (BH or NS) formed by the gravitational 
collapse in Lovelock theory of gravity. In section 2, we present the; bricif overview of generalized Vaidya solution 
in general Lovelock theory of gravity. Next we investigate the gravitational collapse in GB gravity, dimensionally 
continued Lovelock gravity and pure Lovelock gravity in sections 3-5. We have discussed a accretion phenomena 
upon the collapsing object in these gravity theories in section 6. Finally, the paper ends with a short discussions 
in section 7. 
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2 Brief Overview of Generalized Vaidya Solution in Lovelock 
Gravity Theory 



The metric ansatz in (n -|- 2)-dimensional spherically symmetric Vaidya space-time can be written as 



ds^ = -f{v, r)dv^ + 2dvdr + r^dfll 



(3) 
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where r is the radial coordinate (0 < r < oo) and the null coordinate v (—00 < v < 00) stands for advanced 
Eddington time coordinate and f{v, r) = 1 — "^^^-i"^ where m(v, r) gives the gravitational mass inside the 
sphere of radius r and dfl^ is the line element of a n-dimensional unit sphere. The energy-momentum tensor for 
the Vaidya null radiation in the space-time is given by Tab = j-tlah, where la = (1, 0, 0, 0) and /x is the energy 
density of the null radiation. Recently, Cai et al (2008) have calculated the field equations for Vaidya metric 
and their solutions in Lovelock gravity. Here we write briefly the field equations and solutions from their work. 



The Einstein's tensors for Vaidya metric in Lovelock gravity are given by [Cai et al, 2008] 
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The Einstein's field equation is given by 



Gab = SnGTab 



(6) 



(7) 

Since we have = Q^, so from above equation (7) we must have = T^. Now we assume that the 
spherical part of the energy- momentum tensor has the form T? = aT^ = ctT^, where cr is a constant. 

Now from conservation equation VaT^ = 0, we get the following two equations as [Cai et al, 2008] 
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o^t: + drTi + -t; = (8) 
9.t; + ^^iin^T; = (9) 



and 



For pure null radiation, we get = = 0. Here we consider = T^ ^ 0. So from equation (9), we get 
[Cai et al, 2008] 



where C{v) is a function of v. Now define a function 

1 - f{v, r) 
Fiv, r) = 

From (4), (7), (10) and (11), we get the equation 
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which integrates to yield 
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where m{v) is arbitrary function of v, fin = r(i+i^) Q{r) = J r^'^dr i.e.. 
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- — -rr- when a ^ — — 
e(r) = \ (14) 
Inr when cr — —-^ 



Now from equations (5), (7) and (11) , we get 

1^ nl 



Now using (13) and (15), we obtain 

T:^M = #4r + ^^ (16) 

We see that (10) and (16) satisfy (8). From (6), (7) and (11), we get (using T^' = 0) 

z=0 

So the energy-momentum tensor can be written as 



Tab = [llah + {P+P)ilanb + njb) +Pgab (18) 

In the comoving co-ordinates {v,r, 61,62, ...,0ri), the two eigen vectors of energy- momentum tensor namely 
la and Ua are hnearly independent future pointing Ught-Uke vectors (null vectors) having components la = 
(1, 0, 0, 0) and Ua — (//2; —1,0, 0) and they satisfy the relations 

W = nan'' = 0, Ian"" = -1 (19) 

Here, fj, is the energy density of Vaidya null radiation, p is the energy density and p is the radial pressure 
satisfying p = ^ap with p — C{v)r~^^^~"'\ The density p and pressure p come from the Lovelock gravity. The 
solution (13) is called the generalized Vaidya solution in Lovelock theory. Clearly this equation denotes the 
linear barotropic equation of state. We must follow that when cr = — ^ the radiation era is signified whereas 
for (7 = pressureless dust filled model is realised. When a > ^, it denotes dark energy, to be very particular 
quintessence era. At last when cr 1 it implies ACDM whereas a > 1 means phantom era. 



3 Gravitational Collapse in GB Gravity 

Setting p = 2, Co = 0, ci = 1 and C2 = a, the generalized Vaidya solution (13) in Lovelock gravity reduces to 
generalized Vaidya solution in GB gravity. Here the parameter a is called the GB coupling parameter having 
dimension (length)^ and is related to string tension as a^-'^.In the GB theory the function f{v, r) can be 
expressed as [Cai et al, 2008] 
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n — 3 



(20) 



We shall discuss the existence of NS in generalized Vaidya space-time by studying radial null geodesies. In 
fact, we shall examine whether it is possible to have outgoing radial null geodesies which were terminated in 
the past at the central singularity r = 0. The nature of the singularity (NS or BH) can be characterized by the 
existence of radial null geodesies emerging from the singularity. The singularity is at least locally naked if there 
exist such geodesies and if no such geodesies exist it is a BH. 



The equation for outgoing radial null geodesies can be obtained from equation ([3]) by putting ds^ 
dnl = as 



and 
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dv 



(21) 



dr f{v, r)' 

It can be seen easily that r = 0, v = corresponds to a singularity of the above differential equation. 
Suppose X = ^ then we shall study the limiting behaviour of the function X as we approach the singularity at 
r = 0, V = along the radial null geodesic. If we denote the limiting value by Xq then 
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Using (20) and (22), we have 
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Now choosing m{v) = mof " ^ and C{v) = Cqv"'^^ '^^ the equation (23) yields 
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On simplifying wc have the algebraic equation of Xq as 
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- Xo^ + 4Xo - 4 = 



(25) 



This is a very complicated algebraic equation. So as it will be bit difficult to express the value of the root 
Xq in terms of all the parameters, we will vary the value of these parameters to check whether the concerned 
equation gives positive real root for these set of parametric values. Below we have prepared a table 1 where 
some particular sets of values of the parameters have been considered. At first we have to note that as we 
increase the dimension the tendency of having positive root increases. Besides as the value of mo is been 
increased for lower dimension we may not have positive root. So in this case we will have a BH singularity. 
The parameter a takes a leading role in having positive roots. If we take smaller values of a the tendency of 
'not having' positive roots increases even if for larger dimensions too. a plays a role of multiplier to the GB 
term Rgb{= R"'^'^^ Ra^-yS - 'iR°'^Rafi + -R^)) in the action [Boulware et al, 1985]. When a is very low it is 
very obvious that the impact of the GB term in the gravity will be lower. As a — >■ we get back the Einstein 
gravity. So we can predict that tendency of formation of BH reduces due to expanding nature of universe. 
Hence physically we can interpret this phenomenon as : in Einstein gravity circumstances were ready to give a 
BH or to wrap a singularity by as event horizon. On the contrary when we introduce the effect of GB term, i.e., 
we take the expanding universe in account it is observed that the event horizon is not been formed under the 
same circumstances. The negative pressure here is forcing the singularity to be naked. Co has no such impact 
on getting positive roots. 

At last we will discuss the impact of the most important parameter a which is actually determining the 
nature of the fluid present inside the universe. When u > ^ we know it denotes the dark energy (DE) era (both 
quintessence and phantom) as this will not obey the strong energy condition (p + 3p > 0). For this case wc have 
seen that the chance of getting a BH is less in the higher dimensions. In some previous works when Debnath 
et al (2003, 2004), Banerjee et al (2003) have studied the collapse in higher dimensional Tolman-Bondi model 
they have shown that in Einstein gravity the possibility of getting a BH increases in higher dimensions. So we 
have reached an opposite result in our current work i.e., when GB gravity is considered. This may be a strong 
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impact of negative pressure of the fluid. 

So ultimately we may tell that the possibility of the formation of BH in expanding universe is bleak. In lit- 
erature, Babichev et al (2004) and Biswas et al (2011) had shown that the accretion of dark fluid either reduces 
the mass of the BH or weakens its accretion procedure. Similarly, here we can conclude that while forming a 
singularity by collapsing, DE/expanding universe opposes the formation of the BH and rather it increases the 
preferences of growing a NS. 
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Table I: Nature of the roots {Xq) of the equation (25) for various values of parameters involved. 



4 Gravitational Collapse in Dimensionally Continued Lovelock 
Gravity 

If we consider dimensionally continued Lovelock gravity, the solution of the equation (13) yields [Cai et al, 2008; 
Ilha et al, 1997, 1999; Nozawa, 2006] 
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In the first solution of equation put to(w) — mov and C{v) — Cqw""'' and using (22), we have 
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On simplifying we have the algebraic equation of Xq for even values of n as 
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In the second solution of equation (|26p . put m{v) = rngf and C{v) — Cqv '"^ ^ and using (22) and applying 
the limits, we get the algebraic equation of Xq for odd values of n as 
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Like the previous section here also it is tough to determine the explicit form of the solution in terms of 
the different governing parameters. So here also we have varied the values of the parameters and dimension to 
observe their impact on forming positive or non-positive/non-real solution. Like the previous section here also 
we have listed all the data in a table formate (table II) . Now analyzing the table we can say that a has no big 
impact on the nature of root. But if toq and Co are both small then irrespective of the dimension we will get 
positive roots. But as either of these two parameters are increased we will have real positive solution only if 
the dimension is odd. So in this case, the singularity will be naked. For even dimension we will have black holes. 

It is to be noted that when cr — 2, i.e., in phantom era no BH forms. As like the previous section we 
interpret this phenomenon as follows : the negative pressure of DE is opposing a singularity to be a BH and 
forcing it to be a NS. 
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Table II: Nature of the roots (A"o) of the equations (28) and (29) for various values of parameters involved. 
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5 Gravitational Collapse in Pure Lovelock Gravity 

In pure Lovelock gravity, only two coefBcients cq and Cfc are non- vanishing with 1 < fc < [^^] • Choosing 
Cfe = a^*^"^, the solution of f{v, r) can be found from (13) and is given by [Cai et al, 2008] 



{1 - fiv, r)}' = --^^ + ^ I ^ '2:' + ) (30) 
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cor^fc 167rG / m{v) C{v)e{r) 



n(n + l) n vO„r"-2'=+i rn-2fc+i 
where a is a constant length scale. 

Assuming m{v) = mov"~^'^+^ and C{v) = Cqu""^^""'^, using (22) and applying the limits we have, 

1 _ A\^2.-2 ^ levrG jmoTO^jO^r^ ^ g^^n-2.-n. I ^^^^ 



and simplifying, we get the algebraic equation of Xq as 

ISTT^-tG r f 1 + ^) moXo"-2'=+i + IMl^Xo"-^'^-"- - f 1 - na^/^-^ = q (32) 

V 2 / na +1 \ Xq J 

As this is also a complicated equation to give an explicit solution of Xq as a function of the parameters 
and dimension, we will again check the values of roots extracted from the equation for constant values of the 
parameters. All the parametric values, dimensions and the roots found are given below in a tabular form 
(table III). Here we can see unlike the previous cases irrespective of whatever parametric values used we are 
having non-positive roots in each cases. It implies that in pure Lovelock gravity we will always have a black hole. 
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0.11 2(phantom) 0.001 

1 2(phantom) 0.001 

0.3 2(phantom) 0.001 

0.3 .75 (quintessence) 1 
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1 .75 (quintessence) 1 

0.3 0.2 2 

0.11 0.2 2 

1 0.2 2 

1 O(dust) 2 

1 -^(radiation) 2 



Table III: Nature of the roots {Xq) of the equation (32) for various values of parameters involved. 



6 Gravitational collapse when some fluid is also accreting upon the 
collapsing object 

It is an obvious fact that a highly massive star accretes fluid around it mainly from the dust cloud in which 
the star is present or from any companion star which is present in a binary system with the super massive star 
as a companion. As the star absorbs mass the gravitational piill increases irrespective of the increment of the 
electronic force at its surface. So when the system collapses under its own gravity there may be an impact of the 
accreting fluid. Even when the accreting fluid is of DE type this may cause a chance of evolving a NS whereas 
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the fluid like normal matter may increase the chance of forming an event horizon. In this section we will check 
whether our speculation is true or not. 



It is well known that the rate of change of mass of the central object for accreting phenomena is f2„r"TQ^ 
(generalizing the result of Babichev et al (2004)). So the total mass will be 

M - Mdv = M- Qur^'T^dv 

Here, is non-diagonal stress energy tensor component evolved due to the accreting fluid given by (16). So 
equation (20) will be changed to 



fiv, r) = f{v, r) 



fir. 
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+ C(w)e(r) 



(33) 



and equation (22) becomes 



lira 
r 



(34) 



/(^. + ^{^ + c(^)e(r)}] 



where f{v,r) is given in cq.(20). 



For all gravity, we take m{v) = mov"~^^ and C{v) = Cqi^"^^ '^^ we get 

Clearly in the equation (35), Xq is not depending at all upon the GB parameter a or any other information 
of the Vaidya metric given by eq.(20). So if we try to modify eqs.(28), (29) and (32) we will have the same 
equation (35). While the fluid is accreting upon a collapsing object, the nature of singularity (NS or BH) will 
be same for above mentioned three types of gravities. 

Like the previous sections it is really tough to get the explicit solution of Xq as a function of parame- 
ters and dimension. So we will again check the values of roots extracted from the equation for constant values 
of parameters. All the parametric values, dimensions and the roots found are given in a tabular form (table IV). 

Here we have seen for all the ranges of different parameters, irrespective of dimensions we get NS. As we 
have stated before in 2004, Babichev et al have shown that while our universe is expanding the existence of 
highly compact objects like BH at Big Rip is quite impossible. They have shown that the fluid (DE) which is 
responsible for creating such a negative pressure and responsible for the accelerated expansion of present day 
universe decreases the mass of BH while it is been accreted upon the BH (such that before arriving the Big Rip 
point all the BHs will be evaporated completely). Likewise when wc have assumed the fluid is accreted upon a 
collapsing object it will be very obvious for the fluid that it will oppose the collapsing object to have an event 
horizon (i.e., to form a BH) and will have a tendency to produce a NS. 
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mo 


a 


Co 






Positive roots {Xq) 












4D 


5D 


6D 


7D 


8D 


9D 


0.2 


2 (phantom) 


0.5 


1.2455 


1.1555 


1.1034 


1.0688 


1.0431 


1.022 


2 


2 (phantom) 


0.005 


0.6157 


0.645 


0.6729 


0.6969 


0.717 


0.7348 


0.002 


l(ACDM) 


0.005 


1.955 


1.896 


1.807 


1.7051 


1.6090 


1.5270 


0.2 


l(ACDM) 


5 


1.1612 


1.0452 


1.00 


0.9741 


0.9586 


0.9488 


0.2 


0.5(quintessence) 5 


0.2251 


0.2665 


0.3193 


0.3752 


0.4321 


0.4892 


0.02 


0.5 (quintessence 


) 0.05 


1.9408 


1.8700 


1.7730 


1.579 


1.52 


1.5041 


0.01 


0.25 


0.005 


1.7902 


1.6337 


1.4981 


1.398 


1.326 


1.2759 


2.5 


0.1 


0.5 


0.4871 


0.512 


0.5570 


0.5932 


0.6273 


0.6677 



Table IV: Nature of the roots (Xq) of the equation (35) for various values of parameters involved. 



7 Conclusions 

The generalized Vaidya solution in Lovelock theory of gravity in (n + 2)-dimensions have been thoroughly as- 
sumed in this work. Gauss-Bonnet gravity, dimensionally continued Lovelock gravity and pure Lovelock gravity 
have been considered and it has been successfully shown that these three particular forms of Lovelock theory 
of gravity can be constructed by suitable choice of parameters. We have studied the occurrence of singularities 
formed by the gravitational collapse in the above three particular forms of Lovelock theory of gravity. 

In GB gravity three major remarkable results have been observed : (i) Controlling the parameter a when 
we go back to Einstein gravity we find that the possibility of forming a NS increases, (ii) When we increase 
dimensions possibility of finding NS increases while for Einstein gravity contrary to the increase in chances of 
forming NS, the BH is formed for higher dimensions, (iii) The equation of state parameter (a here), when 
denotes negative pressure, tendency of forming NS increases. From the above three observations we can infer 
that for modified gravity theory chances of forming naked singularity is much much higher than forming a 
BH. Dimensionally continued Lovelock gravity theory resembles with the GB gravity in the matter of collapse 
whereas pure Lovelock gravity allows only BH to be formed. We can interpret this as, the gravity in case of 
pure Lovelock gravity theory is more stronger than the expanding force (dark force). At last in section 6 we 
have considered collapse under accretion where we have seen it is almost impossible to form a BH while the 
accretion is going upon the collapsing object in the expanding universe (presented by the modified gravity the- 
ory). Babichcv et al (2004) have shown that under the accretion of phantom fluid the BH mass gets evaporated 
and the concerned BH will never be able to face the Big Rip. But here in the current work we can speculate 
that in expanding universe (specially in phantom era) no BH at all will be formed if we are in Gauss-Bonnet 
gravity or considering accretion procedure upon the collapsing object. 
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